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Abstract 

Polarized microwave emission from dust is an important foreground that may con- 
taminate polarized CMB studies unless carefully accounted for. Modeling of polar- 
ization from dust requires a quantitative understanding of grain alignment. I review 
the current status of grain alignment theory outlining recent advances in quanti- 
tative description of the alignment. In particular, I show that the grain-alignment 
theory is a predictive one, and its results nicely match observations. Those indi- 
cate that the most important process of alignment is related to radiative torques 
acting on irregular grains. The recently developed analytical model of radiative 
torque alignment has proven to be a very efficient tool for predicting the degree of 
grain alignment. We expect the alignment theory to further mature before CMBPol 
flight, which would ensure a better accounting for the dust-related polarization. At 
the same time, CMBPol should provide the additional testing of grain alignment, 
clarifying the reliability of polarimetry for tracing of magnetic field. 



1 Introduction 



Diffuse Galactic microwave emission, apart from carrying important infor- 
mation on the fundamental properties of the interstellar medium, interferes 
with Cosmic Microwave Background (CMB) experiments (see Bouchet et al. 
1999). Polarization of the CMB provides information about the Universe that 
is not contained in the temperature data alone. In particular, it offers a unique 
way to specifically trace the primordial perturbations of tensorial nature (i.e. 
cosmological gravitational waves), and allows one to break some important de- 
generacies that remain in the measurement of cosmological parameters with 
intensity alone (see Zaldarriaga et al. 1997). 

Among different sources of polarized foregrounds, interstellar dust is probably 
the most difficult to deal with (see Lazarian & Finkbeiner 2003 for a review). 
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Polarized radiation from dust arises from grain being aligne Grain align- 
ment has proven to be a very difficult problem. It is well known that the degree 
of grain alignment depends both on grain size and composition, leading to the 
frequency dependence of the polarization (see Kim & Martin 1994, Hildebrand 
et al. 2000). Moreover, in addition to the "classical" power- law distribution 
of larger grains (Mathis, Rumpl & Nordsieck 1977), dust has a population of 
tiny grains (Leger & Puget 1984), which are frequently called PAH. While 
the studies of the alignment of such a particles, which are essentially large 
molecules, is only starting, the theory of alignment of " classical" large grains 
has a long and exciting history (see Lazarian 2003 for a review). 

Theory of grain alignment has been developing fast recently, which makes even 
the most recent review by Lazarian (2007) outdated in a number of respects. 
In what follows we provide an updated outlook on the grain alignment theory 
keeping in mind two interrelated issues, namely, 

a) how grain alignment theory can help in understanding of polarization to be 
measured by CMB-dedicated missions, e.g. CMBPol; 

b) how CMBPol mission can help in understanding of grain alignment. 

There are historic reasons why the researchers may be reluctant to rely on grain 
alignment theory while modeling polarization from dust. Grain alignment has 
a reputation of being an astrophysical problem of the longest standing. In- 
deed, for a long time since the discovery of dust-induced starlight extinction 
polarization in 1949 (Hall 1949; Hiltner 1949) the mechanism of alignment 
was both enigmatic and illusive. Thus, for many years grain alignment the- 
ory did not have predictive powers. Works by great minds like Lyman Spitzer 
and Edward Purcell moved the field forward in terms of understanding the 
basic grain dynamics physics. However, the theory was constantly failing to 
account for new observational data. Fig. la demonstrates the complexity of 
grain motion as we understand it now. 

The weakness of the theory caused a somewhat cynical approach to it among 
some of the polarimetry practitioners who preferred to be guided in their work 
by the following rule of thumb: All grains are always aligned and the align- 
ment happens with the longer grain axes perpendicular to magnetic field. This 
simple recipe was shattered, however, by observational data which indicated 
that 

I. Grains of sizes smaller than a critical size are either not aligned or marginally 
aligned (Mathis 1986, Kim & Martin 1995). 

II. Carbonaceous grains are not aligned, while silicate grains are aligned (see 
Mathis 1986). 



In the absence of the alignment a marginal degree of polarization can arise from 
difference of grain temperatures stemming from the difference in grain cross sections 
in respect to the anisotropic radiation. In what follows, we disregard the latter effect. 



2 



III. A substantial part of small grains grains deep within molecular clouds are 
not aligned (Goodman et al. 1995, Lazarian, Goodman & Myers 1997, Cho & 
Lazarian 2005 and references therein). 

These facts were eloquent enough to persuade even the most sceptical types 
that the interpretation of interstellar polarimetric data does require an ade- 
quate theory. The challenges related to separating of the contribution from 
dust from the one by polarized CMB makes it even more urgent to develop a 
predictive theory of grain alignment. 

One should note that adequate grain alignment theory is able to give boost to 
many astrophysical fields. First of all, polarization from aligned dust is a widely 
used way to study magnetic fields in interstellar medium and molecular clouds. 
Needless to say that ambiguities in grain alignment result in ambiguities of 
the interpretation of polarization in terms of the underlying magnetic fields. 
Moreover, there are numerous astrophysical conditions when dust is bound 
to be aligned, but this fact is under-appreciated. For instance, a common 
explanation of light polarization from comets or circumstellar regions is based 
on light scattering by randomly oriented particles (see Bastien 1988 for a 
review) . The low efficiency and slow rates of alignment were sometimes quoted 
to justify such an approach. However, it has been proven recently that grain 
alignment is an efficient and rapid process. Therefore, we do expect to have 
circumstellar, interplanetary, and cometary dust aligned. Of particular interest 
in this respect are T-Tauri accretion disks (see Cho & Lazarian 2007). Tracing 
magnetic fields in these environments with aligned grains opens new exciting 
avenues for polarimetry. 

In this review I shall show that the modern grain alignment theory is con- 
sistent with the observations available. Moreover, the theory allows us for the 
first time ever make quantitative predictions of the polarization degree from 
various astrophysical objects. However, the latter part of the studies is in its 
infancy. We expect that the growing interest to dust polarization from the 
CMB community will stimulate more work both in grain alignment theory 
and its observational testing. We feel that the CMBpol mission should both 
benefit from our current better understanding of grain alignment and pro- 
vide the critical testing of the remaining controvercies related to the grain 
alignment theory. 

Below, in §2 I shall show how the properties of polarized radiation are related 
to the statistics of aligned grains. In §3 I shall briefly discuss the main align- 
ment mechanisms. In §4 I shall discuss the effects that are related to internal 
dissipation of energy in a wobbling grains. In §5 I shall discuss the predictions 
obtained for the dominant mechanism of alignment, namely, the radiative 
torque alignment. In §6, I shall discuss the prospects of grain alignment the- 
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ory and in §7 I shall discuss the synergy of studying of CMB polarization and 
improving of our knowledge of grain alignment. 



2 Aligned Grains & Polarized Radiation 

A practical interest in aligned grains arises from the fact that their align- 
ment results in polarization of the starlight passing through dust as well as in 
polarization in grain emission. Below we discuss why this happens. 

2.1 Polarized Starlight from Aligned Grains 

For an ensemble of aligned grains the degrees of extinction in the directions 
perpendicular and parallel to the direction of alignment are different!" 2 "! There- 
fore initially unpolarized starlight acquires polarization while passing through 
a volume with aligned grains (see Fig. 2a). If the extinction in the direction 
of alignment is rii and in the perpendicular direction is tj_ one can write the 
polarization, P a b s , by selective extinction of grains as 



where the latter approximation is valid for t\\—t± <C 1. To relate the difference 
of extinctions to the properties of aligned grains one can take into account the 
fact that the extinction is proportional to the product of the grain density 
and their cross sections. If a cloud is composed of identical aligned grains rn 
and Tj_ are proportional to the number of grains along the light path times 
the corresponding cross sections, which are, respectively, C\\ and C±. 

In reality one has to consider additional complications (like, say, incomplete 
grain alignment) and variations in the direction of the alignment axis rel- 
ative to the line of sight. (In most cases the alignment axis coincides with 
the direction of magnetic field.) To obtain an adequate description, one can 
(see Roberge & Lazarian 1999) consider an electromagnetic wave propagating 
along the line of sight (the z° axis, as on Fig. lb). The transfer equations 
for the Stokes parameters depend on the cross sections C xo and C yo for lin- 
early polarized waves with the electric vector, E, along the x° and y° axes 
perpendicular to z° (Lee & Draine 1985). 



2 According to Hildebrand & Dragovan (1995), the best fit of the grain properties 
corresponds to oblate grains with the ratio of axis about 2/3. 
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Fig. 1. (b) Right panel- Geometry of observations (after Roberge & Lazarian 1999). 

To calculate C xo and C yo , one transforms the components of E to the princi- 
pal axes of the grain, and takes the appropriately-weighted sum of the cross 
sections C\\ and C± for E polarized along the grain axes (Fig lb illustrates 
the geometry of observations). When the transformation is carried out and the 
resulting expressions are averaged over the precession angles, one finds (see 
transformations in Lee & Draine 1985 for spheroidal grains, and in Efroimsky 
2002a for the general case) that the mean cross sections are 

C xo = C avg + (C ± - C M ) (l - 3 cos 2 C) , (2) 

C yo = C avg + - R (C± — C||) , (3) 

( being the angle between the polarization axis and the x° y° plane, and 
C avg = (2C j_ + C||) /3 being the effective cross section for randomly-oriented 
grains. To characterize the alignment, we used in eq. (3) the Rayleigh reduction 
factor (Greenberg 1968) defined as 

R= (G(cos 2 6)G(cos 2 (3)) , (4) 

where angular brackets denote ensemble averaging, G(x) = 3/2 (x — 1/3), 9 
is the angle between the axis of the largest moment of inertia (henceforth 
the axis of maximal inertia) and the magnetic field B, while (3 is the angle 
between the angular momentum J and B. To characterize the alignment with 
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respect to the magnetic field, the measures Qx = (G(6)) and Qj = (G(j3)) are 
employed. Unfortunately, these statistics are not independent and therefore R 
is not equal to QjQx (see Lazarian 1998, Roberge & Lazarian 1999). This 
considerably complicates the description of the alignment process. 



2.2 Polarized Emission from Aligned Grains 



Aligned grains emit polarized radiation (see Fig. 2b). The difference in t\\ and 
r± for aligned dust results in the emission polarization: 

p _ (l-e-'i)-(l-e^) ^ n -r ± 
em ~ (1 - e-n) + (1 - e -rx) ~ T|| + Tj _ ' W 



where both optical depths r|| are t± were assumed to be small. Taking into ac- 
count that both P em and P a b s are functions of the wavelength A and combining 
eqs.(l) and (6), one obtains for r = (ry + t±)/2 

P em (A) » -P ate (A)/r(A) , (6) 



which establishes the relation between the polarizations in emission and ab- 
sorption. The minus sign in eq (6) reflects the fact that emission and absorp- 
tion polarizations are orthogonal. This relation enables one to predict the far 
infrared polarization of emitted light if the starlight polarization is measured. 
This opens interesting prospects of predicting the foreground polarization aris- 
ing from emitting dust using the starlight polarization measurements (Cho & 
Lazarian 2002, 2003). As P abs depends on R, P em also depends on the Rayleigh 
reduction factor. 



2.3 Circular Polarization from Aligned Grains 



One way of obtaining circular polarization is to have a magnetic field that 
varies along the line of sight (Martin 1972). Passing through one cloud with 
aligned dust the light becomes partially linearly polarized. On passing the 
second cloud with dust gets aligned in a different direction. Hence, the light 
gets circular polarization. Literature study shows that this effect is well re- 
membered (see Menard et al 1988), while another process entailing circular 
polarization is frequently forgotten. We mean the process of single scattering 
of light on aligned particles. An electromagnetic wave interacting with a sin- 
gle grain coherently excites dipoles parallel and perpendicular to the grain's 
long axis. In the presence of adsorption, these dipoles get phase shift, thus 
giving rise to circular polarization. This polarization can be observed from an 
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Fig. 2. (a)Left panel- Polarization of starlight passing through a cloud of aligned 
dust grains. The direction of polarization (E) is parallel to the plane of the sky 
direction of magnetic field, (b ) Right panel- Polarization of radiation from a optically 
thin cloud of aligned dust grains. The direction of polarization (E) is perpendicular 
to the plane of the sky direction of magnetic field. 

ensemble of grains if these are aligned. The intensity of circularly polarized 
component of radiation emerging via scattering of radiation with k wavenum- 
ber on small (a <C A) spheroidal particles is (Schmidt 1972) 

V(e, e , d) = -^^(anal - ajjaa) ([e x ei]e) (e e) , (7) 



where e and ex are the unit vectors in the directions of incident and scattered 
radiation, e is the direction along the aligned axes of spheroids; a± and an are 
the particle polarizabilities along e and perpendicular to it. 

The intensity of the circularly polarized radiation scattered in the volume 
AT(d, r) at |d| from the star at a distance |r| from the observer is (Dolginov 
& Mytrophanov 1976) 

AV(d,r) = - ^"Z^ , 2 fl([dxr]ft)(dr)Ar(d,r) , (8) 
bi d r d — rr 



where L+ is the stellar luminosity, n^ nst is the number of dust grains per a 
unit volume, and ay is the cross section for producing circular polarization, 
which for small grains is oy = i/ (2k 4 )(a\\a* ± — ajjajj. According to Dolgi- 
nov & Mytrophanov (1976) circular polarization arising from single scattering 
on aligned grains can be as high as several percent for metallic or graphite 
particles, which is much more than one may expect from varying magnetic 
field direction along the line of sight (Martin 1971). In the latter case, the 
linear polarization produced by one layer of aligned grains passes through 
another layer where alignment direction is different. If on passing through a 
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Simplified Model of Alignment 

~b. v Nutation of c5 around J 
in < 10" 5 s 




I J Alignment of co and J 
' co /internal ali gnme nt/ 



~j Precession of J around B 






Gradual alignment of J 



Precession of J is rapid => magnetic field 

is the axis of alignment 
II and ± are possible 




Fig. 3. Left panel- Alignment of grains implies several alignment processes acting 
simultaneously and covering various timescales. Internal alignment was introduced 
by Purcell (1979) and was assumed to be a slow process. Lazarian & Draine (1999a) 
showed that the internal alignment is f 6 times faster if nuclear spins are accounted 
for. The time scale of J and B alignment is given for diffuse interstellar medium. It is 
faster in circumstellar regions and for cometary dust. Right panel- Grain dynamics 
as seen in the grain frame of reference. The Barnett magnetization is directed along 
$7, and it causes a gradual grain remagnetization as the precesses in the grain 



axes. 



single layer, the linear polarization degree is p, then passing through two layers 
produces circular polarization that does not exceed p 2 . 



3 Grain Alignment Theory: Major Mechanisms 



3.1 Characteristic time scales involved in grain alignment 



We have seen in the previous sections that both linear and circular polariza- 
tions depend on the degree of grain alignment given by the Rayleigh reduction 
factor (see Eq. (4)). Therefore it is the goal of the grain alignment theory to 
determine this factor. Table 1 shows that the wide range of different time 
scales involved makes the brute force numerical approach doomed. 

Grain motion as it undergoes alignment is rather complex (see Fig. 3). This 
makes the problem of grain alignment very difficult. 
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A number of different mechanisms that produce grain alignment has been 
developed by now. Dealing with a particular situation one has to identify the 
dominant alignment process. Therefore it is essential to understand different 
mechanisms. 



Table 1 

Time-scales relevant for grain alignment 



Symbol 



Meaning 



Definition 



Value (s) 



tnucl 
tc 
t L 



t-gas 
t E 



VK e h 2 (h-l)J 2 

nuclear relaxation time (^^(j^^Bar 



t ro t thermal rotational period 
tear Barnett relaxation time 



6 x io-st-J/V^V 2 



crossover time 
Larmor precession time 



tRT Radiative precession time 



gas damping time 
electric precession time 



tDG paramagnetic damping time 



27rpc d J|| 
X'VB 

\d<j)/dt\ 

4I \\ 

nmv th b A 

2-k 

n E 

2pa? 
5K(uj)T 2 B 2 



g MxlQ e i ^ r)fl{s y_ 5i JA ) 2 F{Tel 



21.35p 2 a 7 _ 5 / 1 ( S )^/i- 2 (^) 2 F(T ri ) 
L6 X 109 (wft%) 1/2 / 2 ( s )«-5 



GtpT d \„2 2 



1.1 x Kri ^ Jal^s 



1 "W^ 1 



Qe3 ' " ' Au rad ' 
12/ 



4.6 x 10 12 ( T ^f 7 2)s6_ 5 



0.2 x 10 11 p- 1 E- 1 p6ja^ 5 



10 13 B^K^a^s 2 



Notations: 
a: minor axis 
a_5 — a/10~ 5 cm 

h — ratio of moment inertia 

T g — T g /85K : normalized gas temperature 
T ro t = (Tg +T^)/2: rotation temperature 
h — n/20cm ^ : normalized gas density 

X — 10 — ^x/^d'- rea l part of magnetic susceptibility 
K e , n uj: imaginary part of magnetic susceptibility by electron 
fj,^: grain magnetic moment 
Tn — 9n ^" n : magnctogyric ratio nuclei 
p, n — fj, n j : normalized magnetic moment of nucleus 



b: major axis 
s = a/b < 1: ratio of axes 
p — p/3gcm~ 3 : normalized grain density 
= T^/15K: normalized dust temperature 

B — B/5fJ.G: normalized magnetic field 



k e = K e /10~ 



l (T e ) 



Jd 
t- 1 



|| J-L^c 
t\ -I i 



cl 



grain angular momentum at T - 
-f- t 1 ,: total nuclear relaxation time 



cl' 



^rad ~ u rad/ u I 3 RF : energy density of radiation field 

Qe3 — Qr-e3/10 — 2 : third component of radiative torques 

p — 10 15 TJa _5 k € : electric dipole moment 

U — U/0.3V: normalized voltage 

L^ z : magnitude of H2 torque 

W = W/0.2: kinetic energy of H 2 

F(r) = [l + (Q T /2) 2 ] 2 

T e [: electron spin-spin relaxation rate 

h(s) = s- 6 (l + s 2 ) 2 



and nuclear spin 
7e — 9e ^ B '■ magnetogyric ratio for electron 

p, N = eh/2m p c = 5.05 X 10~ 24 ergs G^ 1 

J therm' g ram angular momentum at T — T gas 

can also include inelastic relaxation 

A — A/1.2/j,m: wavelength of radiation field 
E = E/10~ 5 Vcm~ 1 : electric field 
K e — K e /10 2 : electric constant 
Cj — uj / 10° r ad a -1 : angular velocity 
C — C/0.2 fraction of absorbed atoms 
a — a/10 11 cm -2 : density of recombination sites 
T n : nuclear spin-spin relaxation rate 
Me ~ Mb = eh/2m e c: Bohr magneton 
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3.2 Paramagnetic Alignment 



The Davis-Greenstein (1951) mechanism (henceforth D-G mechanism) is based 
on the paramagnetic dissipation that is experienced by a rotating grain. Para- 
magnetic materials contain unpaired electrons that get oriented by the inter- 
stellar magnetic field B. The orientation of spins causes grain magnetization 
and the latter varies as the vector of magnetization rotates in the grain body 
coordinates. This causes paramagnetic loses at the expense of the grain ro- 

— *■ 

tation energy. Be mindful, that if the grain rotational velocity f2 is parallel 
to B, the grain magnetization does not change with time and therefore no 
dissipation takes place. Thus the paramagnetic dissipation acts to decrease 
the component of fi perpendicular to B and one may expect that eventually 
grains will tend to rotate with f2||B provided that the time of relaxation £d_g 
is much shorter than the time of randomization through chaotic gaseous bom- 
bardment, t gas . In practice, the last condition is difficult to satisfy. It is clear 
from Table 1 that for 10~ 5 cm grains in the diffuse interstellar medium, t^-G 
is of the order of 10 13 a^_ 5 ^s 2 B^s , while t gas is 5 x 10 12 n^ 2 o)T^2^ 2 a^ 5 ) s if 
magnetic field is 10~ 5 G and temperature and density of gas are 100 K and 
20 cm -3 , respectively. 

The first detailed analytical treatment of the problem of D-G alignment was 
given by Jones & Spitzer (1967) who described the alignment of J using the 
Fokker-Planck equation. This approach allowed them to account for mag- 
netization fluctuations within the grain material, and thus provided a more 
accurate picture of the J alignment. The first numerical treatment of D-G 
alignment was presented by Purcell (1969). By that time, it became clear 
that the original D-G mechanism is too weak to explain the observed grain 
alignment. However, Jones & Spitzer (1967) noticed that if interstellar grains 
contain superparamagnetic, ferro- or ferrimagnetic inclusions^, the to-G ma y 
be reduced by orders of magnitude. Since 10% of atoms in interstellar dust 
are iron, the formation of magnetic clusters in grains was not far fetched (see 
Martin 1995). However, detailed calculations in Roberge & Lazarian (1999) 
showed that the degree of alignment achievable cannot account for the ob- 
served polarization coming from molecular clouds if grains rotate thermally. 
This is the consequence of the thermal suppression of paramagnetic alignment 
first discussed by Jones & Spitzer (1967). These internal magnetic fluctua- 
tions randomize grain orientation with respect to the magnetic field if the 
grain body temperature is close to the rotational one. 

Purcell (1979) pointed out that fast rotating grains are immune to both 
gaseous and internal magnetic randomization. Substantial degrees of align- 



3 The evidence for such inclusions was found much later through the study of 
interstellar dust particles captured in the atmosphere (Bradley 1994). 
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ment could be possible if the rotation in the same direction could be preserved 
by the external uncompensated torques. Both the radiative torques and other, 
subtle processes that we discuss dealing with grain rotation limit the ability of 
grain alignment in accordance with the Purcell (1979) suggestion (see Spitzer 
& McGlynn 1979, Lazarian & Draine 1999a, Hoang & Lazarian 2008). 

Lazarian & Draine (2000) predicted that PAH-type particles can be aligned 
paramagnetically due to the relaxation that is faster than the DG predictions. 
In fact, they showed that the DG alignment is not applicable to very swiftly 
rotating particles, for which the Barnett magnetic field gets comparable to 
magnetic fields induced by uncompensated spins in the paramagnetic mate- 
rial. For such grains, this relaxation is more efficient than the one considered 
by Davis & Greenstein (1951). This effect, that is termed "resonance relax- 
ation" in Lazarian & Draine (2000), allows the alignment of PAHs. These tiny 
"spinning" grains are responsible for the anomalous foreground microwave 
emission (Draine & Lazarian 1998, see also Lazarian & Finkbeiner 2003 for a 
review) . 



3.3 Mechanical Alignment of Regular Grains: supersonic alignment 

The Gold (1951) mechanism is a process of mechanical alignment of regular 
grains. We shall understand by regular grains the grains that do not exhibit 
helicity (cf. §3.4). Grains presented by spheroids and bricks are the kind of 
regular grains we are talking about. 

The original Gold (1951) idea is easy to understand. Consider a needle-like 
grain interacting with a stream of atoms. Assuming that collisions are inelastic, 
it is easy to see that every bombarding atom deposits with the grain an angular 
momentum 53 = m atom r x v atom , which is directed perpendicular to both the 
needle axis r and the velocity of atoms v atom . It is obvious that the resulting 
grain angular momenta will be in the plane perpendicular to the direction of 
the stream. It is also easy to see that this type of alignment will be efficient 
only if the flow is supersonicF"!. 

The quantitative numerical study of the Gold alignment in Roberge et al. 
(1995) was done under the assumption of the perfect coupling of J with the 
axis of maximal inertia. This study shows a good correspondence with an 
analytical formulae for the alignment of J vector in L94 when the gas-grain 
velocities are hypersonic. An analytical study in Lazarian (1997) accounts for 
the incomplete internal alignment in a more sophisticated way, compared to 



4 Otherwise grains see atoms coming not from one direction, but from a wide cone 
of directions (see Lazarian 1997a) and the efficiency of alignment decreases. 



11 



L94, and predicts the Rayleigh reduction factors of 20% and more for grains 
interacting with the Alfven waves. 

As we discuss below, Purcell (1979) introduced the concept of pinwheel torques 
that act to induce grain rotation with velocities much in excess of the thermal 
Brownian rate. Such fast rotating grains seemed not to be susceptible to the 
Gold (1951) process. This persuaded many researchers that mechanical align- 
ment is of marginal importance. Indeed, it seems natural to accept that fast 
rotation makes it difficult for gaseous bombardment to align grains. However, 
the actual situation is more interesting. First of all, it was proven that mechan- 
ical alignment of suprathermally rotating grains is possible (Lazarian 1995). 
Two mechanisms that were termed "crossover" and "cross section" alignment 
were introduced there. The mechanisms were further elaborated and quantified 
in Lazarian & Efroimsky (1996), Lazarian, Ozik & Efroimsky (1996), Efroim- 
sky (2002b). Second, Lazarian & Drane (1999ab) showed that grains may be 
thermally trapped and not rotate at high speed in spite of the presence of pin- 
wheel torques (see more in §5.3). The original Gold (1951) alignment should 
be applicable to such grains. Moreover, supersonic velocities may be available 
in the interstellar medium due to MHD turbulence (Lazarian & Yan 2002, 
Yan & Lazarian 2003, Yan, Lazarian & Draine 2004). 

Therefore, the problem of the mechanical alignment of regular grains is not 
that it is not efficient, but that in many cases the process may be subdom- 
inant. This is the case, for instance, when a grain has intrinsic helicity. The 
mechanical alignment may be important for small grains, e.g. PAHs, if they 
move supersonically in respect to gas or flow of ions. 

3.4 Mechanical Alignment of Helical Grains: subsonic alignment 

Astrophysical grains are known to be irregular. Therefore they can, for in- 
stance, scatter more impinging atoms to the right than to the left, the same 
way that a model grain in Fig. 4 scatters photons. If present, the helicity 
radically changes the dynamics of grains. 

As far as we know, the mechanical alignment of helical grains was first briefly 
discussed in Lazarian (1995) and Lazarian, Goodman & Myers (1997), but 
was not elaborated there. Lazarian & Hoang (2007b) calculated the torques 
acting on a model helical grains and showed that effective alignment is possible 
for subsonic drifts of grains in respect to the ambient gas. The estimates 
obtained their are very encouraging, but they appeal to the details of grain- 
gas interaction that should be clarified by the further research. 

It is easy to understand why supersonic drift may not be necessary for helical 
grains. For such grains the momentum deposited by regular torques scales in 
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proportion to the number of collisions, while the randomization adds up only 
as a random walk. In fact, the difference between the mechanical alignment of 
spheroidal and helical grains is similar to the difference between the Harwit 
(1971) alignment by stochastic absorption of photons and the radiative torque 
alignment (see § 3.5). Indeed, while the Harwit (1971) alignment requires 
very special conditions to overpower the gaseous and emission randomization 
(Purcell & Spitzer 1971), the radiative torques acting on a helical grain easily 
beat randomization. 

3.5 Radiative Torque Alignment 



Fig. 4. (a) Left panel- A model of a "helical" grain, that consists of a spheroidal 
grain with an inclined mirror attached to it, reproduces well the radiative torques 
(from Lazarian & Hoang 2007a). (b) Right panel - The "scattering coordinate sys- 
tem" which illustrates the definition of torque components: ai is directed along the 
maximal inertia axis of the grain; k is the direction of radiation. The projections 
of normalized radiative torques Q e i, Qe2 and Q e ^ are calculated in this reference 
frame for $ = 0. 

Stochastic radiative torques arising from the absorption from a star were con- 
sidered by Harwit (1971), who addressed the issue of interaction of symmetric, 
e.g. spheroidal, grains with a radiative flow. These torques were shown to be 
subdominant in most astrophysical conditions (Purcell Sz Spitzer 1971). In 
what follows, we do not consider stochastic radiative torques. 

Regular radiative torques were introduced by Dolginov and Mytrophanov 
(1976) who considered "twisted grains" that can be characterized by some he- 
licity. They noticed that "helical" grains would scatter differently the left- and 
right-polarized light, for which reason an ordinary unpolarized light would spin 
them up. The subset of the "helical" grains was not properly identified, and the 
later researchers could assume that it is limited to special shapes/materials. 



Modeling scattering of photon by a grain consisting 
of an ellipsoid and a mirror attached at itstegd. 



View from the right 
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One way or another, this ground-breaking work did not make much impact to 
the field until Draine & Weingartner (1996) numerically showed that grains 
of rather arbitrary irregular shapes experience strong radiative torques. This 
initially was interpreted as the evidence of the fast grain rotation, but later 
research that radiative torques can both spin-up and slow down grains. In 
particular, Lazarian & Hoang (2007a) noticed that grain rotation may be es- 
sentially stopped by radiative torques. In the presence of thermal fluctuations 
grains do not completely stop, but rotate at subthermal rates at low- J attrac- 
tor points (Weingartner & Draine 2003, Hoang & Lazarian 2008a). 

A key point is that anisotropic starlight radiation can align grains on their 
own, i.e. without the influence of any other, e.g. paramagnetic, mechanism. 
In fact, in most cases, in particular in diffuse interstellar gas, the radiative 
torques are so strong that they absolutely dominate the dynamics of grains. 

The effect of alignment induced by radiative torques was discovered by Dol- 
ginov & Mytrophanov (1976). In their paper they considered a tilted oblate 
grain with the helicity axes coinciding with the axis of maximal inertia, as 
well as a tilted prolate grain for which the two axes were perpendicular. They 
concluded, that subjected to a radiation flux, the tilted oblate grain will be 
aligned with longer axes perpendicular to magnetic field, while the tilted pro- 
late grain will be aligned with the longer axes parallel to magnetic field. The 
problem was revisited by Lazarian (1995), who took into account the internal 
relaxation of energy within grains and concluded that both prolate and oblate 
grains will be aligned with longer axes perpendicular to the magnetic field. 
However, Lazarian (1995) did not produce quantitative calculations and un- 
derestimated the relative importance of radiative torque alignment compared 
to other mechanisms. Curiously enough is that our recent analysis of the align- 
ment of an irregular grain mechanical properties of which are presented by a 
triaxial ellipsoid revealed that, even without internal relaxation, grains tend to 
get aligned with long axes perpendicular to magnetic field (Hoang & Lazarian 
2008c). 

The explosion of interest to the radiative torques we owe to Bruce Draine, 
who realized that the torques can be treated with the DDA code by Draine 
& Flatau (1994) and modified the code correspondingly. The magnitude of 
torques were found to be substantial and present for grains of all irregular 
shapes studied in Draine (1996), Draine & Weingartner (1996, 1997). After 
that it became impossible to ignore the radiative torque alignment. More 
recently, radiative torques have been studied in laboratory conditions (Abbas 
et al. 2004). In what follows we identify the alignment by radiative torques with 
the alignment of grains by the anisotropic radiation flux. This corresponds to 
the original thinking in Dolginov & Mytorphanov (1976). 

The Draine & Weingartner (1996, 1997) papers signified a qualitative change 
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in the landscape of grain alignment theory. Although the numerical modeling 
there was rather simplified_J, these papers conjectured that radiative torques 
alignment may be the dominant mechanism in the diffuse interstellar medium. 
Weingartner & Draine (2003) attempted to make more realistic modeling by 
accounting for thermal fluctuations within grains (see §4). 




Fig. 5. Examples of irregular shapes studied in Lazarian & Hoang 2007a. 



DDA, a eff = 0.2 /xm, \=1.2/zm 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

COS© COS© 

Fig. 6. (a) Left panel- Two components of the radiative torques are shown for our 
analytical model (solid lines) in Fig. 7a and for an irregular grain in Fig. 8 (dashed 
lines), (b) Right panel - Radiative torques for different grain shapes. From Lazarian 
& Hoang (2007a). 

The analytical model for radiative torque alignment was proposed in Lazar- 
ian & Hoang (2007a). This simple model was shown to reproduces well the 
essential basic properties of radiative torques. The model consists of an ellip- 
soidal grain with a mirror attached to its side (see Fig. 4). This model allows 
an analytical treatment and provides a physical insight why irregular grains 
get aligned. Fig. 7 illustrates why radiative torques tend to align grains the 
"right way", i.e. in agreement with observations. Interstellar grains experience 
internal relaxation that tends to make them rotate about their axis of max- 
imal moment of inertia. Therefore, it is sufficient to follow the dynamics of 



5 In fact, the crossovers induced by radiative troques were treated incorrectly in 
Draine & Weingartner (1997), which induced a notion of cyclic trajectories. A rig- 
orous treatment of crossovers in Lazarian & Hoang (2007a) corrected this point and 
showed that cyclic trajectories for grains do not exist. Instead, in the absence of 
gaseous bombardment, grains get settled either at an attractor point. 
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angular momentum to determine grain axes alignment. Let us call the compo- 
nent of torque parallel to J the spin-up torque H and perpendicular to J the 
alignment torque F. The angular momentum J is precessing about magnetic 
field due to the magnetic moment of a grain (see Dolginov & Mytrophanov 
1976). The alignment torques F are perpendicular to J and therefore as J gets 
parallel to B the fast precession of the grain makes the torques averaged over 
J precession vanish as £ — > 0. Thus the positions corresponding to J aligned 
with B are stationary points, irrespectively of the functional forms of radiative 
torques, i.e. of components Q e i(Q) and Q e 2(©)- In other words, a grain can 
stay aligned^ with £ = or ir. 

As J, due to internal relaxation (see §4) the grains tend to rotate with J 
aligned with grain shortest axes (corresponding to the maximal moment of 
inertia), the rotation with long grain axes perpendicular to magnetic field 
corresponds to the stationary point. If at the stationary point the spin-up 
torques are positive, the grain rotates fast and we have a high-J attractor 
point. If at the stationary point the the torques spin the grain down, we have 
a point of slow rotation, or low-J attractor point. Interestingly enough, the 
alignment happens also in the absence of internal relaxation (see Hoang & 
Lazarian 2008c). 

Note, that grains can be both "left-handed" and "right-handed". The two 
species of different handedness exhibit their sets of Q e \ and Q e 2 functions 
that are described by the analytical model (see Lazarian & Hoang 2007a). 
For our grain model to become "right handed" the mirror should be turned 
by 90 degrees. The introduction of the analytical model in Lazarian & Hoang 
(2007a) ushered in the quantitative theoretical modeling of radiative torque 
alignment. However, to understand it one should recall the key features of 
grain dynamics related to grain internal relaxation. This is done in the next 
section. The visual comparison of the radiative torques acting on arbitrarily 
shaped grains shown in Fig. 5 and the predictions for the torques arising from 
a model in Fig. 4 is given in Fig. 6. For this comparison we changed the 
handedness of the shape 1 grain (denoted therefore shape 1* in the Fig. 6). 
The phase trajectories of irregular grains and model grains are similar (see 
Fig. 8). 



The arguments above are quite general, but they do not address the question 
whether there are other stationary points, e.g. whether the alignment can also hap- 
pen with J perpendicular to B. To answer this question one should use the actual 
expressions for Q e i(G) and Q e 2{&)- The analysis in Lazarian & Hoang (2007a) 
shows that there is, indeed, a range of angles between the direction of radiation and 
the magnetic field for which grains tend to aligned in a "wrong" way, i.e. with long 
axes parallel to magnetic field. However, this range of angles is rather narrow and 
does not exceed several degrees and, as we discuss further, this range narrow that 
the range of grain thermal wobbling. 
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F is alignment torque ( J_to J) 
H is spin-up torque (// to J) 



Magnetic field 



grain 




Stationary points: 

F=0for 

£ = or 7r 



J/Zfy Fast 



precession 
about B 



F'<0 is an attractor point 
If spin-up torques H>0 we have high-J attractor 



If 



H<0 



low-J attractor 



Fig. 7. A simplified explanation of the grain alignment by radiative torques. The 
grain, which is depicted as a spheroid in the figure, in fact, should be irregular 
to get non-zero radiative torque. The positions J parallel (or anti-parallel) to B 
correspond to the stationary points as at these positions the component of torques 
that changes the alignment angle vanishes. As internal relaxation makes J aligned 
with the axis a\ of the maximal moment of grain inertia, the grain gets aligned with 
long axes perpendicular to B. 



ip = 60 , right helicity, AMO shape t, A = 1 .2/u.m , g = Q.2j&m, tp = 60' 




Fig. 8. Phase trajectory maps are obtained here by assuming that grains have ini- 
tially some large values of angular momentum. The trajectories evolve as a result 
of radiative torque action, (a) Left panel- Phase trajectory map obtained for the 
model grain given shown in Fig. 4. (b) Right panel- The same for an irregular grain 
in Fig. 5 (shape 1). From Lazarian & Hoang (2007a). 



17 



4 Grain dynamics: grains wobbling, flipping and thermally trapped 

4-1 Internal relaxation 

To produce the observed starlight polarization, grains must be aligned, with 
their long axes perpendicular to magnetic field. According to eq. (4) this in- 
volves alignment not only of the grains' angular momenta J with respect to 
the external magnetic field B, but also the alignment of the grains' long axes 
with respect to J (see Fig. 3). The latter alignment is frequently termed the 
"internal alignment". 

Purcell (1979) was the first to consider internal relaxation as the course of 
internal grain axes alignment. He showed that the rates of internal relaxation 
of energy in wobbling grains are much faster than the time of grain alignment. 
This induced many researchers to think that all grains rotate about the axis 
of their maximal inertia and the internal alignment is generically is perfect. 
Lazarian (1994) corrected this conclusion by showing that the fluctuations 
associated with the internal dissipation through the Fluctuation Dissipation 
Theorem should induce the wobbling, which amplitude increases as the effec- 
tive rotational temperature of grains approaches that of grain material. 

The rate of internal relaxation is an important characteristic of grain dynam- 
ics. Spitzer & McGlynn (1979), Lazarian & Draine (1997, 1999ab), Lazarian 
& Hoang (2008), Hoang & Lazarian (2008b) demonstrated its crucial signif- 
icance for various aspects of grain alignment. Purcell (1979) was the first 
to evaluate the rates of internal relaxation, taking into account the inelastic 
effectdj] and a particular effect that he discovered himself, the Barnett relax- 
ation. As we know, the Barnett effect is the magnetization of a paramagnetic 
body as the result of its rotation. Purcell (1979) noticed that the wobbling 
rotating body would experience changes of the magnetization arising from the 
changes of the direction of rotation in grain axes. This would entail fast relax- 
ation, with a characteristic time about a year for a 10~ 5 cm grains. Lazarian 

6 Draine (1999a) identified 10 6 times stronger relaxation related to nuclear 
spins and Lazarian & Hoang (2008) showed that superparamagnetic grains, 
e.g. grains with magnetic inclusions (see Jones & Spitzer 1967, Mathis 1986, 
Martin 1994), exhibit the enhanced internal relaxation. 

The relative role of internal relaxation for the Barnett, nuclear and inelastic 
relaxation is shown in Fig. 9. The calculations of inelastic relaxation from 
Lazarian & Efroimsky (1999) are used for the plot. These rates are important 
for grain dynamics. In particular, they influence how grain crossovers happen 

7 A more recent study of inelastic relaxation is provided in Lazarian & Efroimsky 
(1999). 
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Fig. 9. The relative role of the Barnett, nuclear and inelastic relaxation for grains 
of different sizes. Rates of Barnett relaxation are from Purcell (1979), nuclear relax- 
ation are from Lazarian & Draine (1999b), inelastic relaxation are from Lazarian & 
Efroimsky (1999). 

(see §5.3). 



4-2 Regular pinwheel torques 



Purcell (1975, 1979) realized that grains may rotate at a rate much faster 
that the thermal rate if they are subject to systematic torques. P79 identified 
three separate systematic torque mechanisms: inelastic scattering of impinging 
atoms when gas and grain temperatures differ, photoelectric emission, and H2 
formation on grain surfaces (see Fig. 10). Below we shall refer to the latter as 
"Purcell's torques". These were shown to dominate the other two for typical 
conditions in the diffuse ISM (P79). The existence of systematic H 2 torques 
is expected due to the random distribution over the grain surface of catalytic 
sites of H 2 formation, since each active site acts as a minute thruster emitting 
newly-formed H2 molecules. A later study of uncompensated torques in Hoang 
& Lazarian (2008b) added additional systematic torques to the list, namely, 
torques arising from plasma-grain interactions and torques arising from emis- 
sion of radiation by an irregular grains. Radiative torques arising from the 
interaction of the isotropic radiation with an irregular grain (Draine & Wein- 
gartner 1996) also represent the systematic torques fixed in the grain body. 
We shall call all the systematic torques above pinwheel torques to distinguish 
them from the systematic torques arising from anisotropic flow of photons or 
atoms interacting with an helical grain (see below). 
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Spin-down 



Torques decrease J! 



Spin-up 
Torques 
increase J|| 



Fig. 10. (a) Left panel- A grain acted upon by Purcell's torques before and after a 
flipover event. As the grain flips, the direction of torques alters. The H2 formation 
sites act as thrusters. (b) Right panel- A regular crossover event as described by 
Spitzer <fe McGlynn (1979). The systematic torques nullify the amplitude of the 
J component parallel to the axis of maximal inertia, while preserving the other 
component, J±. If J± is small then the grain is susceptible to randomization during 
crossovers. The direction of J is preserved in the absence of random bombardment. 



Purcell (1979) considered changes of the grain surface properties and noted 
that those should stochastically change the direction (in body-coordinates) of 
the systematic torques. Spitzer & McGlynn (1979, henceforth SM79) devel- 
oped a theory of such crossovers. During a crossover, the grain slows down, 
flips, and thereafter is accelerated again (see Fig. 10). 

From the viewpoint of the grain-alignment theory, the important question is 
whether or not a grain gets randomized during a crossover. If the value of the 
angular momentum is small during the crossover, the grains are susceptible 
to randomization arising from atomic bombardment. The original calculations 
in Spitzer & McGlynn (1979) obtained only marginal correlation between the 
values of the angular momentum before and after a crossover, but their anal- 
ysis disregarded thermal fluctuations within the grain with temperature Td us t- 
According to the Fluctuation-Dissipation Theorem these thermal fluctuation 
induce thermal wobbling with the frequency determined by the internal relax- 
ation rate (see Lazarian 1994). If the crossover happens over time larger than 
the grain internal relaxation time, Lazarian & Draine (1997) showed that the 
Spitzer & McGlynn (1979) theory of crossovers should be modified to include 
the value of thermal angular momentum J d ^ (2kT dust I\\) 1 ^ 2 , where iy is the 
maximal moment of inertial of an oblate grain. What is the size of grains for 
which the effect gets important? When nuclear relaxation is accounted for, 
this provides a size of a c ~ 10~ 4 cm. Thus for large grains, e.g. in dark clouds 
and accretion disks the modification of the Spitzer & McGlynn (1979) theory 
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0.1 1.0 10.0 

°cri( 1 0~ 5 cm) 

Fig. 11. The applicability range for the Spitzer & McGlynn (1979), Lazarian <fe 
Draine (1997) and Lazarian & Draine (1999ab) theories. The latter predicts flipping 
the rate of which were elaborated in Hoang & Lazarian (2008b). J m ax parameterizes 
here the value of the pinwheel torques, i.e. J max = t gas dJ/dt, where the derivative 
dJ/dt arises from the pinwheel torque action. 

is important. 

4-3 Thermal trapping of grains 

What would happen for small grains? Lazarian & Draine (1999a) predicted 
that such grains would flip fast and this will alter the Purcell torques making 
grain rotate thermally in spite of the presence of the uncompensated torques 
(see Fig. 10). This way Lazarian & Draine (1999a) explained why interstellar 
grains smaller than 10 -6 cm poorly aligned via paramagnetic relaxation. A 
more elaborate study of the flipping phenomenon in Roberge & Ford (2000 
preprint; see also Roberge 2004) supported this conclusion. Weingartner (2008) 
challenged this result by showing that spontaneous thermal nipping does not 
happen if the internal relaxation diffusion coefficient (Lazarian & Roberge 
1997) is obtained with a different integration constant. Does this invalidate 
the phenomenon of thermal trapping predicted in Lazarian & Draine (1999a)? 
Our study in Hoang & Lazarian (2008b) shows that the flipping and trapping 
do happen, if one takes into account additional fluctuations associated, for 
instance, with the action of the uncompensated torques. 

What does happen with the grains which do not flip, but smaller than a c ? 
For those grains the original SM79 theory of crossovers is applicable and these 
grains are marginally aligned via Davis- Greenstein process, provided that the 



21 



pinwheel torques are short-lived, i.e. the time scale of their existence is much 
smaller than the time of paramagnetic relaxation. Such grains rotate at high 
rate in accordance with the Purcell (1979) model. As we discuss further, the 
fact that pinwheel torques are not suppressed by flipping may allow better 
alignment by radiative torques. 

Fig. 11 defines the range over which the different model of crossovers are appli- 
cable. It is clear that for a sufficiently large J ma x the flipping gets suppressed. 



5 Quantitative modeling of radiative torques 



5. 1 Predictions of the analytical model 



How can one characterize radiative torques acting on grains of arbitrary shapes? 
Our work showed that the entire description of alignment may be obtained 
with the two components of the radiative torques Q e \ and Q e2 as they are 
defined in the caption of Fig. 4. The third component Q e s is responsible for 
grain precession only. The functional dependences of the torque components 
that are experienced by our model grain are similar to those experienced by 
irregular grains shown in Fig. 5. It is really remarkable that our model and 
grains of very different shapes have very similar functional dependences of 
their torque components (see Fig. 6)! 

For typical conditions of the diffuse ISM the radiative torques are strong and 
absolutely dominate the grain dynamics. In this situation, the alignment de- 
pends only on the ratio of the two components of torques, namely, on Q e i/Qe 2 - 
This ratio is a function of the angle between the radiation direction, but ,as 
the functional dependences of Q e \ and Q e2 for irregular grains correspond well 
to the analytical expressions obtained for a helical model grain in one can 
clearly see from Fig 13, where the function 

(A 2 )(Q e2 ) = / [QT 9Ular (&) - Q? 2 odel (&)] 2 dQ (9) 



characterizes the deviation of the calculated torques Q e2 calculated numeri- 
cally for irregular grains from the analytical prediction of the torques in Lazar- 
ian & Hoang (2007a) model. 

It is very important that one can fully characterize the ratio of these two 
torques by defining the ratio of their amplitude values. In other words, the 
knowledge of the ratio q = Q™f x /Q™ 2 ax i s sufficient to characterize the align- 
ment by radiative torques. This is the second remarkable fact on the properties 
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Fig. 12. The magnitude of the ratio q = Q™^^ /Q^ ax that characterize the radiative 
torque alignment of grains depends on both grain shape and the wavelength of 
radiation. From Lazarian & Hoang 2007a. 



of radiative torques. In terms of practical calculations, this enormously simpli- 
fies the calculations of radiative torques: instead of calculating two functions 
Qei(©) and Q e 2(©) it is enough to calculate just two values Q™£" x and Q^ ax . 
According to Lazarian & Hoang (2007a) the maximal value of the function 
Q e i(Q) is achieved for = of the function Q e 2(0) is achieved at O = ir/4. 
In other words, one can use a single ratio q = Q™^ /Q 1 ^ = Q e i(0)/Q e 2( 7r /4) 
instead of two functions to characterize grain alignment. Thus, it is possible 
to claim that the g-ratio is as important for the alignment as the grain axis 
ratio is important for producing polarized radiation by aligned grains. 

What does affect the g-ratio? One can see that it is not universal. The geo- 
metric optics calculations for the model grain in Fig. 4 show that this ratio 
varies with the change of the angle that characterizes the twist of the mirror in 
respect to the axis of inertia of the grai Similarly, we can see in Fig. 12 that 
the g-ratio is a function of both the grain shape and radiation wavelength. A 
general tendency is that q peaks when grains are of the size comparable to the 
radiation wavelength and decreases below unity when grains are substantially 



8 This change somewhat alters the functional dependences of the radiative torques 
and we find that the best agreement can be obtained by fixing the angle of the 
mirror twist at 7r/4. 
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Fig. 13. Numerical comparison of the torques calculated with DDSCAT for irregular 
grains for different wavelength and the analytical model (AMO) of a helical grain 
(from Lazarian & Hoang 2007a). 

smaller than the radiation that they interact. 



The practical use of the g-ratio is presented in Fig. 16 where the parameter 
space for the existence of only low-J as well as coexistence of both high-J 
and low-J attractor points is shown. As we discussed earlier, in the case when 
high-J attractor point exists, grains eventually get into that position, which 
is characterized by both the fast rotation and perfect alignment in respect to 
magnetic field. In the situation when only low-J attractor point exists, the 
alignment in respect to magnetic field is not stable and the characteristic de- 
grees of alignment vary on average from 20% to 50%. Comparing results in 
Fig. 12 and 16 one may conclude that in the situations when the alignment 
arises from the radiation with the wavelength A larger than the grain char- 
acteristic size, the perfect alignment is likely when the angle ip between the 
radiation direction and magnetic field is larger than n/ 4. 

Are other alignment directions feasible? Due to the fast Larmor precession of 
grains, the magnetic field constitute for them the symmetry axis and therefore 
the alignment can happen either with long axes perpendicular to magnetic 
field ("right alignment) or with long axes parallel to magnetic field ("wrong 
alignment"). The model in Lazarian & Hoang (2007a) shows that for a range 
of angles between the radiation and the magnetic field the alignment gets 
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Fig. 14. Parameter space for which grains have only low- J attractor point and both 
low- J and high attractor point (from Lazarian & Hoang 2007a). In the situation 
when the high- J attractor point is present grains eventually get there and demon- 
strate perfect alignment. In the situation when only low- J attractor point is present, 
the alignment is partial. 

"wrong" , i.e. with the long axes parallel to magnetic field. However, this range 
is narrow (limited to radiation direction nearly perpendicular to magnetic 
field). Therefore the range of "wrong" alignment disappears in the presence 
of grain thermal wobbling described in §4 as was shown in Hoang & Lazarian 
(2008a). 

All the calculations above assumed that the radiation comes from a single 
direction. Is it possible to provide prescriptions for the alignment when the 
radiation field is characterized by smooth distribution of light coming from 
different directions. In this case it is natural to decompose the radiation field 
into multipoles and consider the alignment by individual components. This 
was done in Hoang & Lazarian (2008b) who considered the radiative torque 
alignment by the dipole and quadropole components of the radiation field. 



5.2 Radiative torque alignment of superparamagnetic grains 



Superparamagnetic grains, i.e. grains with enhanced paramagnetic relaxation, 
were invoked by Jones & Spitzer (1967) within the model of paramagnetic 
alignment. What does happen when the dynamics of grains is determined 
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Fig. 15. A paramagnetic grain gets only a low-J attractor point. For the same set 
of parameters a superparamagnetic grain gets also a high-J attractor point. The 
fact that most of the phase trajectories go in the direction of the low-J attractor 
point illustrates the dominance of the radiative torques for the alignment even in 
the case of the superparamagnetic grain. However, high-J attractor points are more 
stable than the low-J attractor points. As a result, all grains eventually end up at 
the high-J attractor point. 

by radiative torques? We see from Fig. 16 that for a substantial part of the 
parameter space grains are driven to the low-J states, i.e. subthermally, which 
is in contrast to the assumption in Draine & Weingartner (1996) that in the 
presence of radiative torques most of the interstellar grains should rotate at 

T ro t ^ Tg as . 

The picture above, however, is different if grains paramagnetically dissipate 
the energy on the time scales shorter than the gaseous damping time. Lazarian 
& Hoang (2008) found that if grains have superparamagnetic inclusions, they 
always get high-J attractor point. Fig. 15 shows that for superparamagnetic 
grains subject to a diffuse interstellar radiation field most grains still get to 
the low attractor point, which reflects the fact that it is the radiative torques 
rather than paramagnetic ones that dominate the alignment. As the high-J 
attractor point is more stable compared to the low-J attractor point, similar 
to the ordinary paramagnetic grains, superparamagnetic grains get transfered 
by gaseous collisions from the low-J to high-J attractor points. Thus, super- 
paramagnetic grains always rotate at high rate in the presence of radiative 
torques. One concludes that, rather unexpectedly, intensive paramagnetic re- 
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Fig. 16. Left panel. Grain alignment by radiative torques in the absence of pinwheel 
torques. Right panel. Grain alignment by radiative torques in the presence of pin- 
wheel torques. The shown case corresponds to the presence of both the low- J and 
high-J attractor points in the absence of pinwheel torques. In the case when only 
a low-J attractor point exists the strong pinwheel torques lift the low-J attractor 
point enhancing the alignment. 

laxation changes the rotational state of the grains, enabling them to rotate 
rapidly. The alignment of grains at high-J point is perfect. 



5.3 Radiative torque alignment in the presence of pinwheel torques 

Pinwheel torques were considered by Purcell (1979) in the context of para- 
magnetic alignment. How do these torques also affect the radiative torque 
alignment? Hoang & Lazarian (2008b) showed that the sufficiently strong 
pinwheel torques can create new high-J attractor points. Therefore for suffi- 
ciently strong pinwheel torques, e.g. for torques arising from H2 formation, one 
may observe the correlation of higher degree of polarization with the atomic 
hydrogen content in the media, provided that H 2 torques as strong as they con- 
sidered in Purcell (1979) and the subsequent papers (see Spitzer & McGlynn 
1979, Lazarian 1995, Lazarian & Draine 1997). Another implicit assumption 
for observing this correlation is that the grains are not superparamagnetic. 
For superparamagnetic grains the alignment, as we discussed above, is perfect 
anyhow. 

The existing uncertainty is related to the magnitude of the pinwheel torques. If 
the correlation is observed, this would mean that, first of all, grains are not su- 
perparamagnetic and, in addition, the pinwheel torques are sufficiently strong 
to overwhelm the radiative torques which drive the grains to low attractor 
points. 
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6 Prospects of Grain Alignment Theory 

6.1 Alignment of grains > 3 x 1CT 6 cm 

Grain alignment theory so far was mostly aimed at understanding of the align- 
ment of large grains, e.g. larger than 3 x 10~ 6 cm, for which we have the 
observational evidence of the alignment. Thus we, first of all, deal with such 
grains. 

Grain alignment theory has witnessed a substantial progress in the recent 
years. For the first time in the history of the subject we can talk about the 
qualitative correspondence of the theoretical predictions and observations (see 
Whittet et al. 2008). Other attempts of polarization modeling were made while 
the quantitative understanding of radiative torques was only at its initial stage, 
but they also provide a satisfactory correspondence with observations (see Cho 
& Lazarian 2005, Pelkonen et al. 2007, Bethell et al. 2007). All this testify 
in favor of radiative torques as a dominant mechanism for the alignment of 
astrophysical grains. 

The questions that require further studies include resolving the issue of whether 
interstellar grains are superparamagnetic or paramagnetic. As we discussed 
above, if grains are superparamagnetic, their alignment is much easier to pre- 
dict. In most of the interstellar conditions we expect such grains to be per- 
fectly aligned with long axes perpendicular to the magnetic field. If, however, 
grains are paramagnetic, they can be sometimes perfectly, if their g-ratio and 
the angle of illumination corresponds to the part of the parameter space in 
Fig. 16 that exhibit a high- J attractor point. Otherwise, a reduced degree of 
alignment is expected. 

How can one distinguish between the two cases above? It looks that two ap- 
proaches are readily available. One is based on statistical studies using a sim- 
ulated distribution of grains with different q parameters and comparing the 
observations and simulations. The other approach is based on identifying of 
the situations when the geometry of illumination is known, as may be the case 
for some circumstellar regions, and studying the alignment there. 

Alternatively, if the correlation of the degree of alignment with the pinwheel 
torques is detected for diffuse interstellar medium (see §5.3) this would also 
be the evidence that the grains are paramagnetic. Indeed, the interstellar 
radiation field is sufficiently strong render strong radiative torques essentially 
everywhere in the diffuse galactic gas. Therefore, the superparamagnetic grains 
are expected to be aligned perfectly according to §5.2 The variations, on the 
contrary, mean that the grain alignment is incomplete. 
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It is also important to identify the possible niches for other alignment pro- 
cesses, in particular to the mechanical alignment of helical grains. As we dis- 
cussed earlier, the model grain in Fig. 4 is helical not only in respect to ra- 
diation, but also to mechanical flows (see also LH07). It is easy to see that 
the functional dependence of the torques does not depend on whether photons 
or atoms are reflected from the mirror. Therefore our extended discussion of 
radiative torque alignment above is applicable to the mechanical alignment of 
helical grains. In the presence of dynamically important magnetic field, the 
alignment is expected with long axes perpendicular to B. If atoms stick to the 
grain surface and then are ejected from the place of their impact, this changes 
the values of torques by a factor of order unity. The helium atoms are likely 
to have specular impacts, however. 

What would it take to make a grain helical for mechanical interactions? This 
is a question similar to one that worried researchers with the radiative torques 
before Bruce Draine modified his code to treat radiative torques (see Draine 
1996). We do not have the simulations of mechanical torques on irregular 
grains, but in analogy with the radiative torques, one could claim that such 
torques should be generic for an irregular grain. Detailed calculations should 
clarify the actual degree of helicity of irregular grains and the magnitude of 
the mechanical torques on such grains. The corresponding calculations are 
straightforward and we expect to have the answer before the CMBPol launch. 



6.2 Alignment of Small Grains 

Radiative torques predict the alignment of interstellar grains larger than ps 
3 x 1CT 6 cm, which is in an excellent agreement with the empirical Serkowski 
(1973) curve (see Lazarian 2007). For particles much less than the wavelength 
the efficiency of radiative torques drops as (a/A) 4 (see Lazarian 1995). Within 
circumstellar regions, where UV flux is enhanced smaller grains can be aligned 
by radiative torques. This could present a possible solution for the reported 
anomalies of polarization in the 2175 A extinction feature (see Anderson et 
al 1996) which have been interpreted as evidence of graphite grain alignment 
(Wolff et al 1997). If this alignment happens in the vicinity of particular 
stars with enhanced UV flux and having graphite grains in their circumstellar 
regions, this may explain why no similar effect is observed along other lines of 
sight. 

The maximum entropy inversion technique in Kim & Martin (1995) indicates 
that grains larger than a particular critical size are aligned. This is consistent 
with our earlier discussion of radiative torques and the Serkowski law (see 
§6.2). However, an interesting feature of the inversion is that it is suggestive 
of smaller grains being partially aligned. Initially, this effect was attributed to 
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the problems with the assumed dielectric constants employed in the inversion, 
but a further analysis that we undertook with Peter Martin indicated that 
the alignment of small grains is real. Indeed, paramagnetic (DG) alignment 
must act on the small grains^. An important feature of this weak alignment 
is that it is proportional to the energy density of magnetic field. This opens 
a way for a new type of magnetic field diagnostics. As very small grains may 
emit polarized radiation as they rotate (see §4.2) both UV and microwave 
polarimetry may be used to estimate the intensities of magnetic field. 



6.3 Alignment of PAHs 

In the range of frequencies from 30GHz to 90GHz the polarized emission from 
PAHs may play an important role. There have been no systematic work on the 
PAHs alignment however. The fact that no polarization from infrared PAH 
emission has been detected does not mean that the PAHs are not aligned. 
Indeed, the infrared emission arises from PAH absorption of a UV photon. 
The likely effect of such an absorption is the randomization of grain axes 
in terms of the angular momentum J. This effect arises from the coupling of 
rotational and vibrational degrees of freedom via internal relaxation (Lazarian 
1994, Lazarian & Roberge 1997, Weingartner 2008, Hoang & Lazarian 2008b). 
Therefore while J can be aligned with magnetic field, the grain axes may be 
randomized in respect to J during the short periods of PAH emission. On 
the contrary, the dipole emission arising from PAHs in accordance with the 
"spinning dust" process in Draine & Lazarian (1998) proceeds all the time 
and, moreover, is relatively insensitive to the randomization of grain axes. 

The only quantitative study of PAH alignment is one in Lazarian & Draine 
(2000). There the so-called "resonance paramagnetic relaxation" process was 
appealed to. The traditional Davis-Greenstein (1951) relaxation is not applica- 
ble to small grains, as it disregards the important magnetization arising from 
the Barnett effect. This magnetization is very strong for the rapidly rotating 
PAHs. Moreover, the magnetization that arises is exactly of the right mag- 
nitude to induce paramagnetic resonance and therefore the dissipation much 
stronger than in the Davis-Greenstein process. Nevertheless, the polarization 
arising from the PAHs predicted in Lazarian & Draine (2000) is marginal for 
the emission at frequencies higher than 30GHz (see Fig. 17) 



To avoid a confusion we should specify that we are talking about grains of 
10~ 6 cm. For those grains the results of Davis-Greenstein relaxation coincide with 
those through resonance relation in Lazarian & Draine (2000). It is for grains of the 
size less than 10 -7 cm that the resonance relaxation is dominant. 
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Fig. 17. Polarization for both resonance relaxation and Davis-Greenstein relaxation 
for grains in the cold interstellar medium as a function of frequency. From Lazarian 
k Draine 2000. 

7 Grain alignment theory and CMBPol 

7. 1 Zodiacal dust as a polarized foreground 

The contribution from the Zodiacal dust is taken into account for the emission 
studies of CMB. At the same time, radiative torque mechanism should make 
the Zodiacal dust at least partially aligned. Indeed, our study in Lazarian & 
Hoang (2007a) shows that the radiative torques do not decrease when grain 
sizes are larger than the wavelength of impinging radiation, contrary to the 
predictions in Dolginov & Mytrophanov (1976). Moreover, Hoang & Lazarian 
(2008c) showed that grains should be aligned even when they are so large 
enough that the internal relaxation becomes negligible. The polarized emission 
by Zodiacal light may influence the studies by CMBPol and therefore should 
be accounted for. 

A possible study of the polarization arising from the Zodiacal light may be 
attempted in far infrared prior to the CMBPol mission. However, even optical 
studies may clarify this issue, if we take into account that scattering from the 
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Fig. 18. Polarization from magnetic grains. Upper panel: Polarization of thermal 
emission from perfectly aligned single domain grains of metallic Fe (solid lines) 
or hypothetical magnetic material that can account for the Foreground X (broken 
lines). Lower panel: Polarization from perfectly aligned grains with Fe inclusions 
(filling factor is 0.03). Grains are ellipsoidal and the result are shown for various 
axial ratios. From Draine & Lazarian 1999. 

aligned grains can result in circularly polarized radiation, as it discussed in §8 



7.2 "Magnetic" dust and "spinning dust" as polarized foregrounds 



The problem discussed in §5.2 was the superparamagnetic nature of the dust. 
Superparamagnetic and, in general, strongly magnetic dust produces polar- 
ized microwave magneto-dipole emission, which can interfere with the CMB 
polarization studies (Draine & Lazarian 1999). If the magnetic properties of 
dust are better constrained as a result of testing of grain alignment predic- 
tions, this should simplify the removal of dust magnetic contribution at the 
microwave range. Similarly, a better understanding of PAHs alignment should 
facilitate the removal of the "spinning dust" polarized contribution. 

The predictions for the polarization arising from magnetic grains are shown 
in Fig. 18 
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7.3 Synergy of CMBPol and grain alignment studies 

Understanding of grain alignment is essential for the successful separating of 
the CMB polarization from the polarization arising from dust. Any construc- 
tion of a polarization template for the dust foreground uses the assumptions 
on grain alignment. 

The progress in understanding of grain alignment can be gauged by model- 
ing of polarization in interstellar gas, molecular clouds, accretion disks. The 
correspondence of the theory and observations obtained along these directions 
is encouraging (see Lazarian 2007). However, more detailed modeling is re- 
quired. We hope that by the time of CMBPol mission the grain alignment 
theory will be better tested and therefore reliable. However, additional tests 
that CMBPol data can provide will also be useful. For the purposes of testing 
grain alignment the larger coverage of frequencies is advantageous. Thus the 
project EPIC-CS or EPIC-2m would be preferable. Quantitative predictions 
of grain alignment that are possible these days (see §5.1) allow quantitative 
testing. 

In other words, CMBPol requires better understanding of grain alignment to 
get cleaner polarization signal from the early Universe and grain alignment 
theory can benefit from the CMBPol mission. The importance of the latter 
should not be underestimated. Indeed, grain alignment theory is aimed at stud- 
ies of magnetic fields in various environments. In Lazarian (2007), apart from 
interstellar medium, other environments, e.g. dark clouds, accretion disks, 
comets, galactic nuclei, AGNs and Seyfet galaxies, are discussed. Therefore 
the progress in understanding how well magnetic fields are traced by aligned 
dust particles is essential for many astrophysical studies. 

Combining polarimetric data at very different wavelength it is possible to 
get insight into magnetic field structure. Examples of such studies combining 
far infrared, near infrared and optical polarimetry are discussed in Lazarian 
(2007). In addition, one may mention yet another interesting possibility. The 
studies of molecular cloud column densities with the near infrared scattered 
light were presented in Padoan et al. (2006) and Juvela et al. (2006). Those 
have shown that large scale mapping of scattered intensity is possible up to 
A v ~ 10 even for clouds illuminated by the average interstellar radiation field. 
The polarization of scattered light should be affected by grain alignment. 
This opens interesting prospects of detailed mapping of magnetic fields at 
sub-arcsecond resolution. 
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